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983704 

The present invention relates to a method for locating a cased bore- 
hole froa an adjacent open or uncased borehole. Many times during the drill* 
ing of oil wells it is necessary or desirable to locate a previously drilled 
well, which is cased with a ferromagnetic material, from a second well which 
is uncased. For example, in the case of a blowout in a previously drilled 
deep well, it may be more economical to intercept the previous well at a 
shallow depth above the total depth and kill the formation than to drill to 
the total depth of the original well. Likewise* when drilling multiple wells 
from a single offshore platform, it is desirable to know the exact location 
10 of the existing wells with respect to the well that is being drilled to pre- 
vent an accidental interception of previously drilled wells. 

In the past it has been impossible to accurately predict the loca- 
tion of a previously drilled well from a second well. While the need for such 
a device has long been felt in the petroleum industry, it has not been satis- 
fied by available equipment. 

Cfcject of the invention is a method for determining the distance 
and direction to a cased borehole from measurements made in adjacent uncased 
borehole, said cased borehole being cased with a material having a remanent 
magnetization. 

20 According to the invention there is provided a method for deter- 

mining the distance and direction to a cased borehole from measurements made 
in adjacent uncased borehole, said cased borehole being cased with a material 
having a remanent magnetization, said method comprising: measuring, in the 
uncased borehole at a plurality of vertically spaced positions, the total 
strength of the magnetic field existing at each of the positions; measuring, 
in the uncased borehole at substantially the same positions, the deflection 
of a magnetic compass from a geographical reference; measuring the depth of 
the spaced positions from a known reference point; assuming a value and 
direction for the earths magnetic field at the positions; and determining 

30 by machine calculations the distance and direction to the cased borehole 

using said measurements and assumed value for the earth's field and location 
for the cased borehole, said calculations using a least squares fit of the 
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measured values to calculated values in the formula for the magnetic induc- 
tion surrounding a long cylinder. 

The present invention provides a method for locating a previously 
drilled well which is cased with a material having a remanent magnetization, 
to have a remanent magnetization in the casing is not difficult since most 
well casing is electromagnetically inspected before it is installed. The 
electromagnetic inspection leaves a remanent magnetization in the casing. 

Since casing is normally installed in individual sections that are 
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Joined together, the remanent magnetization of unperturbed casing normally is 
periodic. 

Magnetometers are avialable for measuring the total strength of a 
magnetic field. The field measured by the magnetometer will, of course, he 
the combination of the field from the magnetised casing plus the field of the 
earth. Knowing the total magnetic field and the earth 1 s field, one can cal- 
culate possible locations of the previously drilled well by assuming various 
locations and computing the minimum least squares error of fit to observed 
total field data. After calculating possible locations of the previously cas- 
10 ed well and knowing the strength of the earth's field, one can determine the 
strength of the field contributed by the cased well. Knowing the location 
and strength of the field contributed by the casing, one can calculate the 
deflection of a magnetic compass caused by the cased well. The calculated 
deflection can then be compared to the actual deflection of the magnetic com- 
pass as determined by a survey that includes both a gyro and magnetic compass. 
Prom a comparison of the calculated deflection with the actual deflection, one 
can then determine the location of the cased well. In those cases where one 
is only interested in knowing the location of the first well so that the second 
well does not intercept it, one can utilise the Information to control the 
20 direction of the second well. 

The present invention will be deacribea by way of example in more 
detail with reference to the attached drawings in which: 

Figure 1 in an isometric view of the locations of the two wells; 
Figure 2 is a representative plot of the total magnetic fiela measur- 
ed by a magnetometer, the magnetic and gyro compass headings} and 

Figures >5 are representative plots of the deflection of a magnetic 
compass ana calculated deflections based on the measured magnetic field. 

The theory of magnetic fields including both the earth's field and 
an additional field is well understood by those skilled in physics. This theory 
50 can be applied to the case of a hollow, infinitely long, magnetise* cylinder 
extending into the earth. When the magnetised cylinder is a well casing, one 
can use the theory of magnetic fields to derive the following formulas for the 
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magnetic field due to a magnetized well casing: 

(1) B n = 5 I & h (2rnip/IQ cos (2nvZ/L) ♦ b h (2rnip/lQ sin (2n*Z/L) 

p n=i n K o (2nWL) n K o (2mia/i} 

<g=*o r r 

(2) B 7 a„ *b(2nirp/L3 sin (2nuZ/L) - b o(2mrp/IQ cos (2mr2/L 

n=! n K o (2nua/L) !T(2nWE) 

In the formula, B p is the radial component of the external field 
and Bg is axial component , a is the outer radius of the casing and L is the 
wave length of the periodic magnetization of the casing. K Q and are modi- 
fied Bessel functions of the second kind, a and b are constants, p is the 

n n 

distance from the axis of the casing to the point of measurement, and Z is 
the vertical distance. 

Formulas 1 and 2 above apply to the case where the magnetized casing 
has a regular periodicity in the magnetization.. Normally, the casing will be 
subjected to only a single stress since the casing is suspended from a hanger 
member in the borehole and remains in tension. In the case of a blowout or 
other abnormality, the casing can be subject to factors that distort the per- 
iodic character of the magnetic field external to the casing. In this case, 
L is chosen larger than the interval over which data is being fit and form- 
ulas 1 and 2 can still be used. 

In Figure 1, there is shown a schematic three-dimensional view of 
the two wells, one having a magnetized liner 10 and the other well 11 being 
drilled to intercept the first. The liner 10 is ferromagnetic and has a re* 
aanent magnetization. Hie remanent magnetization may result from an electro- 
magnetic inspection of the casing before it was installed or may be purposely 
provided by other means. The horizontal coordinate system assumes that the 
first well is located at the axis of the coordinate system and the second well 
is located to the east of the first well and drilled at an angle to inter- 
cept the first well. At known depths in the second well, it is assumed that 
the first well can be in the various positions along a line 12, such as 13, 
while the second well is in the position 18. Thus the distances from the 
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assumed positions of the first well to the second well are H and D or they 
can be represented by the single lino p which is the shortest distance from 
the axis of the second well to the assuaed position of the first well in the 
horizontal plane. The north-south and east-west magnetic axes are also 
shown. 

The method of the present invention consists of the steps of first 
measuring in the second well the total magnetic field. As pointed out above, 
this measurement can best be made using a sensitive magnetometer; for example, 
a magnetometer having a sensitivity of at least l(f 4 gauss would be suitable. 
At the same time that the magnetic field is measured, the particular depths in 
the second well should also be measured. In addition, the deflection of a 
magnetic compass in relation to a gyro compass or similar device should also 
be measured. 

From the total measured field, one can proceed to calculate the 
strength of the radial field B p and vertical field B z resulting from the mag- 
netized casing of the first well at the position of the second well. The 
earth's field is assumed to have a value of about .55 gauss and a direction of 
approximately 65° from the horizontal and directed along an azimuth corres- 
ponding to magnetic north. The exact values for the earth's field will vary 
at various geographical locations. Those given above are approximate values 
for the central Mississippi area of the United States. The total measured 
field can be expressed vectorially as follows: 

wherein B is the total field, B H the earth's horizontal field, B the magnet- 
ized casing's radial field; B y the earth's vertical field; B z the magnetized 
casing's vertical field; a the angle between magnetic North and a line join- 
ing the two wells; and e 1# e 2 and e 3 are unit vectors in the North, East and 
vertically downward directions, respectively. 

If both sides of the above vector equation are squared, one obtains: 



. Cos a + 2B B ♦ 
P v z 



D 2 .B 2 
P z 
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Proa an assumed position for the first well, values for p and a can 

bo obtained and the unknowns c and b in the formulae 1 and 2 calculated. 

n n 

This calculation is done with L greater than the length of the path over which 
the calculations are being considered. Also n should be varied to provide 
sufficient wave lengths to determine the best fit of the observed data to the 
calculated data. The above calculation is carried out by first considering 
only the first order terms, 2 B^B^ Cos a ♦ 2 B y B^ to obtain preliminary value 
of the unknowns a n and b^. These values are then used to calculate the values 
of the squared terms [b h 2 ♦ B y 2 ] ♦ |j> p 2 * B^ 2 ]' From the calculated 
values of the squared terms, new calculations are made for the first order 
terms. The calculations are continued until the best fit of the calculated 
value of the total field and the observed value of the total field is obtain- 
ed. The best fit can be obtained using a conventional least squares fit of 
the neasured data to calculated values, i.e., B 2 measured -B 2 calculated. 

While no methods have been described in detail for carrying out the 
above calculations, anyone skilled in the computer art can program a computer 
to provide a least squares fit of measured data do calculated data. The use 
of least squares to fit measured data to calculated data is described in many 
books on applied mathematics. Thus, no description of a program for a parti- 
cular computer is included. 

After the strengths of the radial and vertical fields in the second 
well resulting from the magnetised liner in the first well have been calcul- 
ated, the resulting deflection of a magnetic compass can be calculated. The 
calculated values of deflection can then be compared to the measured values 
of deflection. This is important since positions of the first well on 
opposite sides of the second well may give substantially the same errors. 
While the least mean square errors are the same, the deflection of magnetic 
compass will be materially different at the two positions. Shown in Figure 3 
is the apparent magnetic declination derived from measurements of the azimuth 
of the second well with magnetic and gyromultishot survey instruments. The 
data are representative of that obtained in a well . Figures 4 and 5 represent 
magnetic declinations calculated for two assumed locations of the first well. 
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It can be seen that the assumed location for tho calculated declinations 
shown in Figure 4 is the location of the first well. 

The above procedure was used to intercept a well which had been 
drilled to approximately 20,000 feet and cased with a material having 
remanent magnetization. The well subsequently blew out and it was desired 
to intersect the well at approximately 10,000 feet to avoid the drilling of 
a relief well to 20,000 feet to kill the formation. Also, it was believed 
that the casing in the original well had parted at a few thousand feet and 
placed the casing at the depth of interest in compression. The data shown 
in Figure 2 represent the nagnetometer data obtained over an interval of 
approximately S00 feet in the second well using a magnetometer of high 
sensitivity. 
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THE EMBODIMENTS OF THE INVENTION IN WHICH AN EXCLUSIVE 
PROPERTY OR PRIVILEGE IS CLMMED ARE DEFINED AS FOLLOWS: 



1. A nethod for determining the distance and direction to a cased bore- 
hole from measurements made in adjacent uncased borehole, said cased borehole 
being cased with a material having a remanent magnetisation, said method com- 
prising: measuring, in the uncased borehole at a plurality of vertically 
spaced positions, the total strength of the magnetic field existing at each 
of the positions; measuring, in the uncased borehole at substantially the 
same positions, the deflection of a magnetic compass from a geographical re- 
ference; measuring the depth of the spaced positions from a known reference 
point; assuming a value and direction from the earths magnetic field at the 
positions; and determining by machine calculations the distance and direction 
to the cased borehole using said measurements and assumed value for the 
earths field and location for the cased borehole, said calculations using a 
least squares fit of the measured values to calculated values in the formula 
for the magnetic induction surrounding a long cylinder, 

2. The method of claim 1 wherein the radial component of the field 
contributed by the casing is given by the formula: 



B . > a *i(2nWL) K. (2mrp/L) 

p n - 1 n K 6 CWL) C0S {2mrZ/L) + b n K o (2nira/L) si * < 2nirZ/L > 

and the axial component of the field contributed by the casing is given by 
the formula: 

B z , a n *o C2ntrp/L 3 sin (2nirZ/L) - b ' 0 »"» /W cos (2n*Z/L) 

n»l n K o (2nira/L) n K q {Wl) 

where a is the outer radius of the casing, Kg and ^ modified Bessel func- 
tions of the second kind, a n and b R constants, I the wave length of the 
periodic magnetisation, p the distance from the axis of the cased hole to 
the point of measurement and Z the distance parallel to the axis. 

3. The method of claim 2 wherein L is the average length of one joint 

of casing. 
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4. The method of claim 2 wherein L is larger than the interval over 
which data is being fit. 

5. The method of claim 2 wherein the measured data is substituted in 
the formula and the least squares fit of the data calculated to the measured 
data obtained, the calculated results then being used to determine the 
calculated deflection of a magnetic compass, 

6. The method of claim 5 wherein the actual deflection of the magnetic 
compass is compared with the calculated deflection to determine the location 
of the cased borehole. 
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